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Abstract 

The optical excitations in Ceo and higher fullerenes, including isomers of C76, C78, and Cs4, are 
theoretically investigated. We use a tight binding model with long-range Coulomb interactions, treated 
by the Hartree-Fock and configuration-interaction methods. We find that the optical excitations in 
the energy region smaller than about 4 eV have most of their amplitudes at the pentagons. The 
oscillator strengths of projected absorption almost accord with those of the total absorption. Next, 
off-resonant third order susceptibilities are investigated. We find that third order susceptibilities of 
higher fullerenes are a few times larger than those of Cgo- The magnitude of nonlinearity increases as 
the optical gap decreases in higher fullerenes. The nonlinearity is nearly proportional to the fourth 
power of the carbon number when the onsite Coulomb repulsion is 2t or At, t being the nearest 
neighbor hopping integral. This result, indicating important roles of Coulomb interactions, agrees 
with quantum chemical calculations of higher fullerenes. 
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^ ■ 1. Phason Lines and Linear Absorption 

Recently, the fullerene family Cat with hollow cage structures has been intensively investigated. A lot 
of optical experiments have been performed, and excitation properties due to 7r-electrons delocalized 
on molecular surfaces have been measured. For example, the optical absorption spectra of Cgo and C70 
Q\ . [1,2] have been reported, and the large optical nonlinearity of Ceo [3,4] has been found. The absorption 
spectra of higher fullerenes (C76, C78, Cs4, etc.) have also been obtained [5,6]. For theoretical studies, 
we have applied a tight binding model [7] to Ceo, and have analyzed the nonlinear optical properties. 
Coulomb interaction effects on the absorption spectra and the optical nonlinearity have been also 
studied [8]. We have found that the linear absorption spectra of Cqo and C70 are well explained by 
the Frenkel exciton picture [9] except for the charge transfer exciton feature around the excitation 
energy 2.8 eV of the Ceo solids [2]. Coulomb interaction effects reduce the magnitude of the optical 
nonlinearity from that of the free electron calculation [8], and thus the intermolecular interaction 
effects have turned out to be important. 

In the previous paper [10], we have extended the calculation of Ceo [9] to one of the higher fullerenes 
C76. We have discussed variations of the optical spectral shape in relation to the symmetry reduction 
from Ceo and C70 to C76: the optical gap decreases and the spectra exhibit a larger number of small 
structures in the dependences on the excitation energy. These properties seem to be natural when we 
take into account of the complex surface patterns composed of pentagons and hexagons. In order to 
understand the patterns clearly, the idea of the phason lines (Fig. 1) has been introduced [11] using 
the projection method on the honeycomb lattice plane [12]. There are twelve pentagons in C76. Six 
of them cluster on the honeycomb lattice, with one hexagon between the neighboring two pentagons. 
There are two groups of the clustered pentagons. The phason line runs as if it divides the two groups. 
We use the following Hamiltonian: 

H = H + H int . (1) 

The first term of eq. (1) is the tight binding model: 

ffo = -* £ (ct CT c i)CT + h.c.), (2) 
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where t is the hopping integral and q j(7 is an annihilation operator of a 7r-electron with spin a at the 
ith carbon atom of the fullerene. It is assumed that t does not depend on the bond length, because 
main contributions come from excitonic effects due to the strong Coulomb potential. The results do 
not change so largely if we consider changes of hopping integrals by bond distortions. We assume the 
following form of Coulomb interactions among 7r-electrons: 



2'* hi •>+ 2' 

+ E W{rij)(£, c l c ^ ~ X )(E C W - 1), (3) 
where r^j is the distance between the ith and jth atoms and 

W{r) = + (r/roVy (4) 

is the Ohno potential. The quantity U is the strength of the on-site interaction, V means the strength 
of the long-range Coulomb interaction, and r is the average bond length. 

The model is treated by the Hartree-Fock approximation and the single excitation configuration 
interaction method, as was used in the previous papers [9,10]. In ref. 9, we have varied the parameters 
of the Coulomb interactions, and have searched for the data which reproduce overall features of 
experiments of Cqq and C70 in solutions. We have found that the common parameters, U = 4t and 
V = 2t, are reasonable. Thus, we use the same parameter set for higher fullerenes. The quantity t is 
about 2 eV as shown in ref. 9. The Coulomb interaction strengths depend on the carbon positions. We 
use the lattice coordinates which are obtained by the public program FULLER [13,14]. The optical 
spectra become anisotropic with respect to the orientation of the molecule against the electric field 
of light, as reported in the free electron model (Hiickel theory) [15]. We obtain numerical optical 
absorption spectra by summing the data of three cases, where the electric field of light is along the 
x-, y-, and z-axes. 

We use a projection operator in order to extract contributions to the optical spectra from a certain 
part of fullerenes. If we write the a projection operator to a part of lattice site set as P, the oscillator 
strength between the ground state \g) and the excited state \k) is written: 

/«,* = E k [\(k\PxP\ 9 )\ 2 

+ \(k\(1 - P)x(l - P)\g)\ 2 

+ (g\PxP\K){K\(l-P)x(l-P)\g) 

+ ( 9 \(1-P)x(1-P)\k){k\PxP\ 9 )}, (5) 

where E K is the excitation energy, and the electric field is parallel with the x-axis. In eq. (5), the first 
term is the contribution from the projected part, and the three other terms are the remaining part. 
The total optical absorption is calculated by the formula: 

+ + ( 6 ) 

K 

where p(uj) = ^/[tt(uj 2 + 7 2 )] is the Lorentzian distribution of the width 7. The projected absorption 
is calculated by eqs. (5) and (6). The projected part does not satisfy a sum rule. So, this results in a 
singularity where excitation energy is large. We will discuss the optical spectra in the energy region 
far from the singularity. 

Figure 2 shows the molecular structures and optical spectra of the C70 molecule and C76 with the 
D2 symmetry, which have been found in experiments. The black atoms are the carbons along the 
phason lines. The hatched circles are the pentagonal carbons. In C70, the phason line runs along 
the ten carbons which are arrayed like a belt around the molecule. In C76, the phason line is located 
almost along the outer edge of the molecule (Fig. 1). The total optical absorption is shown by the 
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bold line, and the absorption from all the pentagonal carbons is shown by the thin line. We find that 
the optical excitations in the energy region lower than It are almost composed of the excitations at the 
pentagonal sites. This property is common to C70 and D2-C76, and also to the T^-G^ for which the 
calculated data are not shown. In higher energy regions, the thin lines give relatively larger oscillator 
strengths, but this is an artifact of the projected wavefunctions. The absorption spectra calculated 
from the projected wavefunctions do not satisfy the sum rule, i.e., the area between the abscissa and 
the curve does not become constant regardless of the excitation wavefunctions. The similar artifact 
will be found in the figures shown afterwards, believe that the projected optical absorption spectra 
are reliable in low energy regions only. Therefore, we limit our comparison of the spectra to the energy 
region lower than about 2t ~ 4 eV. 

In Ceo, the edges of the pentagons are the long bonds, and the bonds between the neighboring 
hexagons are short bonds. The wavefunctions of the fivefold degenerate highest-occupied-molecular- 
orbital (HOMO) have the bonding property, and that the threefold degenerate lowest-unoccupied- 
molecular-orbital (LUMO) has the antibonding property. As the carbon number increases, hexagons 
are inserted among pentagons. The wavefunctions near the LUMO of the higher fullerenes still have 
the antibonding properties, thus they tend to have large amplitudes along the edges of pentagons which 
have the characters like long bonds of Cqq. Recently, the bunching of the six energy levels higher than 
the LUMO has been discussed in the extracted higher fullerenes [16]. The wavefunctions near the 
LUMO distribute on the pentagons. This fact can be understood as the properties characteristic 
to antibonding orbitals. As the excited electron mainly distributes at the pentagonal carbons, the 
electron-hole excitation has large amplitudes at these pentagons. Thus, the oscillator strengths of the 
low energy excitations are mainly determined by wavefunctions at the pentagonal carbons. This is the 
reason why the projected absorptions nearly accord with the total absorptions in the energy regions 
smaller than about 2t. 

If the projections are performed onto each pentagon, we can know contributions to optical spectra 
from the projected carbon sites. We would like to look at this feature, for example, in D2-C7Q. There 
are three carbon atoms, which are not equivalent with respect to symmetries, in this isomer. These 
pentagons are indicated by the symbols, A-C, in Fig. 3(a). The projected absorption spectra are shown 
by thin curves, superposed with the total absorption in Figs. 3 (b-d). The projected absorption is 
multiplied by the factor 12, in order to compare with the total absorption. We find that the projected 
absorption exhibits small structures in the energy region smaller than 2t. The structures depend on 
the kind of carbons. The spectral shapes and oscillator strengths are much far from those of the total 
absorption. It would be difficult to assign experimental features of the total absorption with a set of 
the limited number of carbon atoms. The excitation wavefunctions at the twelve pentagons give rise 
to the shape of the absorption spectra totally. 

2. Nonlinear Optical Response 

In this section, we investigate nonlinear optical properties of higher fullerenes. We focus on the 
off-resonant third order susceptibility in order to estimate the magnitudes of the nonlinear optical 
responses of each isomer. The Coulomb interaction strengths are also changed in a reasonable range, 
because realistic strengths are not well known in higher fullerenes. Based on our results for the optical 
properties of Ceo and C70 [8,9], we can assume V = U/2. The onsite Coulomb strength is varied 
within the range < U < 4t, t being the hopping integral between nearest neighbor carbon atoms. 

In Fig. 4, the relations between the absolute value of the off-resonant susceptibility and the energy 
gap are shown for three Coulomb interaction strengths: U = Ot, 2t, and At. Here, the energy gap is 
defined as the optical excitation energy of the lowest dipole allowed state, in other words, the optical 
gap. For each Coulomb interaction, the plots (squares, circles, or triangles) cluster in a bunch. When 
the energy gap becomes larger, the susceptibility tends to decrease. However, the correlation between 
the susceptibility and the energy gap is far from that of a smooth function. The correlation is merely 
a kind of tendency. Therefore, the decrease in the energy gap of higher fullerenes is one origin of the 
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larger optical nonlinearities of the systems. The actual magnitudes of nonlinearities would also be 
influenced by the detailed electronic structures of isomers. 

In the calculations for Ceo reported previously, the magnitudes of the THG at the energy zero 
are approximately 1 x 1CT 12 esu in the free electron model, and approximately 2 x 1CT 13 esu for 
U = At and V = 2t. In the present calculations for higher fullerenes, the magnitudes are a few 
times larger than those of C%q. Thus, the author predicts that nonlinear optical responses in higher 
fullerenes are generally larger than in Ceo- In our previous paper [8], we discussed the fact that the 
local field correction factor is of the order of 10 for Ceo solids. Since the distance between the surfaces 
of neighboring fullerene molecules in C70 and C76 solids is nearly the same as in Cqq solids, we expect 
that local field enhancement in thin films of higher fullerenes is of a magnitude similar to that in Cgo 
systems. 

TABLE I. Coulomb interaction dependence of the power a where |x (3) (°)l ~ A-N a . 



u 


a 


ot 


5.253 


2t 


4.133 


At 


3.536 



It is of some interests to look at carbon number dependence of the magnitude of the optical 
nonlinearity of the calculated isomers in higher fullerenes. Figure 5 shows |x^(0)| as functions of 
the carbon number N for three Coulomb interaction strengths, U = Ot, 2t, and At. The solid lines 
indicate the linear fitting in the logarithmic scale: |x^ 3 Ho)l ~ A • N a . The powers a for the three 
Coulomb interaction strengths are summarized in TABLE I. When U = Ot, the power a is about 5. 
As U increases, a decreases. It is among 4, when U ~ 2t and At. This magnitude of the power 4 agrees 
with the result of the quantum chemical calculation of higher fullerenes upto Cs4 [17]. Therefore, we 
have shown important roles of Coulomb interactions in nonlinear optical response of higher fullerenes. 

Experimental measurements of optical nonlinearities in higher fullerenes whose carbon number is 
larger than 70 have not been reported so much, possibly because of the difficulty in obtaining samples 
with good quality and the difficult measurements. However, the recent report of the degenerate four- 
wave-mixing measurement of C90 in solutions [18] indicates the larger optical nonlinearity than that in 
Ceo- The magnitude of is about eight times larger than in Ceo, and is apparently enhanced from 
that of the theoretical predictions: (90/60) 4 = 1.5 4 = 5.063. Therefore, further experimental as well 
as theoretical investigations of nonlinear optical properties in higher fullerenes should be fascinating 
among scientists and technologists of the field of photophysics. 
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Figures should be requested to harigaya@etl.go.jp. 
Caption: 

Fig. 1. Phason lines in higher fullerenes. 

Fig. 2. Molecular structures and theoretical optical spectra of (a) D 5 d-C7o and (b) D2-C7Q. In 
the molecules, the black atoms are along the phason lines, and the hatched atoms are the pentagonal 
carbons. In the absorption spectra, the bold line is the total absorption, and the thin line is the 
absorption by the wavefunctions projected on the twelve pentagons. The units of the abscissa are 
taken as arbitrary, and the energy is scaled by t. The parameters are U = 4t, V = 2t, and 7 = 0.06i. 

Fig. 3. (a) The molecular structure of D2-C7Q. The symbols, A-C, indicate the symmetry nonequiv- 
alent pentagons. The figures, (a), (b), and (c), compare the absorption projected on one of the three 
pentagons with the total absorption. The bold line is the total absorption, and the thin line is the 
projected absorption. The units of the abscissa are taken as arbitrary, and the energy is scaled by t. 
The data of the thin line are multiplied by the factor 12. The parameters are U = At, V = 2t, and 
7 = 0.061 

Fig. 4. The absolute value of the off-resonant susceptibility |x^(0)l f° r ^60 an d seven isomers 
of higher fullerenes, plotted against the energy gap (shown in units of t). The squares, circles, and 
triangles represent results for U = Ot, 2t, and At, 

Fig. 5. The absolute value of the off-resonant susceptibility |x^(0)l f° r ^60 and seven isomers of 
higher fullerenes, plotted against the carbon number N. The squares, circles, and triangles represent 
results for U = Ot, 2t, and At, respectively. The left and bottom axes are in the logarithmic scale. The 
solid lines are the results of the linear fitting in the logarithmic scale: \x^\ ~ A ■ N a . 
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